Vortex Bubble Formation in Pair Plasmas 
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It is shown that delocahzed vortex sohtons in relativistic pair plasmas with small temperature 
asymmetries can be unstable for intermediate intensities of the background electromagnetic field. 
Instability leads to the generation of ever-expanding cavitating bubbles in which the electromagnetic 
fields are zero. The existence of such electromagnetic bubbles is demonstrated by qualitative argu- 
ments based on a hydrodynamic analogy, and by numerical solutions of the appropriate Nonlinear 
Schrodinger equation with a saturating nonlinearity. 
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Nonlinear dynamics of vortex solitons, structures that 
permeate different fields of physics, like Nonlinear Optics, 
Bose-Einstein Condensate and Plasma Physics, has re- 
ceived considerable attention in the recent past [T] . The 
vortex solitons, usually, emerge as solutions of a general- 
ized (2+l)D-dimensional nonlinear Schrodinger equation 
(NSE) with a local nonlinearity. The NSE with a de- 
focusing (for instance, cubic) nonlinearity admits stable 
vortex soliton solutions. These solutions, with an angular 
27r phase ramp, and appearing as local dark minima in an 
otherwise bright background, are viewed as the most fun- 
damental two-dimensional solitons embodied in the NSE. 
The NSE with a focusing nonlinearity also supports local- 
ized vortex solitons with phase dislocations surrounded 
by one or many bright rings. The latter, however, are 
unstable against symmetry-breaking perturbations that 
can cause the breakup of rings into filaments. 

New and interesting behavior arises when the non- 
linearity can change its character (as a function of the 
amplitude, for instance) from focusing to defocusing or 
vice versa. A representative example is a sign -changing, 
cubic-quintic (CQ) nonlinearity (focusing for lower and 
defocusing for higher intensities). The NSE with a CQ 
nonlinearity can, under certain conditions, support ex- 
istence of stable localized as well as delocahzed vortex 
soliton solutions (see for details [5]). Though CQ anzats 
is frequently invoked to model light beam propagation in 
various optical media ^4, (like non-Kcrr crystals, chalco- 
genide glasses, semiconductors and so on), its validity to 
describe these materials does require additional explo- 
rations. 

In contrast to the optical media, an NSE with a 
focusing-defocusing nonlineartity can naturally emerge 
from the basic physics models when one investigates the 
propagation of large amplitude electromagnetic pulses in 
a plasma [S]. In some of our recent publications [B],[7], 
we showed that the dynamics of the short intense electro- 
magnetic pulses propagating in a relativistic pair plasmas 
with small temperature asymmetry does, indeed, obey an 



NSE with, what may be called, a saturation nonlinearity. 
Though always positive (unlike the CQ nonlinearity), the 
saturation nonlinearity vanishes for high intensities, and 
thus combines the focusing-defocusing natures. 

The particular form of the saturation function in the 
NSE derived in |S],[7], leads to the existence of localized 
as well as delocahzed optical vortex soliton solution (we 
term these solitons as LOVS and OVS respectively). The 
stability of LOVS (both single and multi-charged) has 
been investigated in detail and the stability areas were 
established. Although the stability of the OVS has been 
established for high intensities of the background fields 
[7] , for intermediate intensities the vortex solitons can be 
unstable. Indeed Kim et al 8 (see also |9j) demonstrated 
numerically that delocahzed vortices described by NSE 
with CQ saturation turns out to be unstable for certain 
intermediate intensities leading to the formation of an 
ever-expanding circular ring. One is naturally faced with 
the question if similar behavior would pertain for the sat- 
uration nonlinearity. Unfortunately the standard linear 
modal stability analysis is inadequate to deal with non- 
exponential growth that could be associated with non- 
Hermitian (non-self-adjoint) linear operators. To deter- 
mine the stability of the OVS in a pair plasma, therefore, 
we resort to a numerical solution of the NSE. 

The nonlinear (2+l)D evolution of an electromagnetic 
(EM) pulse propagating in an arbitrary pair plasma with 
temperature asymmetry can be described by an NSE. 
The dimensionless evolution equation for the slowly vary- 
ing envelope of the vector potential (in co-moving with 
group velocity frame) reads [7]: 
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has a very unusual feature - in the ultra-relativistic Umit 
(|j4p ^ 1) it tends to be and the nonlinear term 
vanishes. 



Before proceeding further, we summarize the steps 
that led to (IT]). We had started with a circularly po- 
larized EM pulse ^ A (ic + iy) cxp(ikoz — cjot) with 
mean frequency ujq and mean wave number fco and 
where A is the slowly varying amplitude. In Eq.(llJ), 
V^ = d^/dx^ + d^/dy^ is the diffraction operator, 
and various dimensionless quantities have been defined 
as: A = |e|l/(mG(To-)c2), r = {ue/cek^^'^) r, t = 

(a;g/2 Wo e^fc)i; where uJe — (47re^no/"i)^/^ is the elec- 
tron Langmuir frequency and m is the electron mass. 
The charges q^ and masses m^ of positive and nega- 
tive ions are assumed to be same (we mainly concentrate 
on the specific case of an electron-positron plasma, i.e. 
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equilibrium state of the system is characterized by an 
overall charge neutrality rig = n^ = no where Uq and 
Uq are the unperturbed number densities of the posi- 
tive and negative species respectively. The background 
temperatures of plasma species are Tq (Tg*" ^ Tq ) and 
mG{z'^) = mK3{z^)/K2{z'^) is the "effective mass", 
[z^ = mc^/T^], where K^ are the modified Bessel func- 
tions. For the nonrelativistic temperatures (T* ^ mc^) 
G* = 1 + 5T*/2toc^ and for the ultra-relativistic tem- 
peratures (r± > TOaC^) G=^ = 4r±/TOc2 > 1. The pa- 
rameter e = [G(ro+)-G(ro-)]/G(ro+) measures the back- 
ground temperature asymmetry of plasma species. For 
the nonrelativistic temperatures e = 5(rQ^ — Tq")/2?7ic^ 
while in ultrarelativistic case e — {Tf^ — Tq)/T^. The 
numerical factor k — 1/4 for non-relativistic tempera- 
tures (= 1/3 for ultrarelativistic temperatures). In de- 
riving Eq.Q with (pi), we have assumed that the plasma 
is highly transparent and the longitudinal extent of the 
pulse is much shorter than its transverse dimensions. The 
medium is self-focusing {df/d\A\ > 0) for \A\ < 1, while 
it becomes defocusing {df/d\A\ < 0) for higher intensi- 
ties (1^1 > 1). 

To investigate the Vortex soliton solutions of Eq. (IT]) , 
one makes the ansatz A = ^(r)exp(iAt -|- imO), where 
the integer to defines the topological charge, and A 
is the nonlinear frequency shift. In polar coordinates, 
this substitution converts Eq. (IT|) into an ordinary dif- 
ferential equation. For a nonzero topological charge, the 
solution has a node at the origin: A^^q — )■ rl^l^o 
where Aq is a constant that measures the slope of A 
at the origin. Depending on the values of A and ^o the 
soliton solution is either localized Ar^oo ^^ or delo- 
calized A^^oo — ^ ^oo (where A = /(^oo) )• The de- 
localized vortex solitons are found to exist only in the 
range < A < 0.25 (oo > A^o > 1) for any value of the 
topological charge m. Thus, while near the vortex core 
the medium is focusing, the background intensity of the 
soliton Aao > 1 i.e., far beyond of the vortex core the 



medium is defocusing and consequently the background 
field is modulationally stable. 

In what follows we consider the stability of OVS. It 
is well known that single charge (to = 1) OVSs are sta- 
ble in self-defocusing media. The multicharged vortices, 
however, despite their very long life [lOj . should even- 
tually decay into m = I vortices [1]. The last state- 
ment (which is based on general topological considera- 
tion) should be valid for any kind of media and will not 
be addressed here. For the particular system described 
by Eq. (IT]), OVSs are stable against azimuthal pertur- 
bation for any values of the allowed background field [7] ; 
the radial stability, on the other hand, can be established 
just for high intensities. Based on the results obtained 
in 8] one should expect that for the intermediate inten- 
sities of the EM field OVSs could be unstable for our 
mixed-type nonlinearity. 

To augment our expectation it is useful to employ 
the hydrodynamic analogy based on so-called Madelung 
transformation A = p^l"^ exp(j0), where p — \A\^ and </> is 
the real phase [TP. With simple algebra (and changing a 
time coordinate as i — > i/2 ), we can transform the NSE 
(fll) into the following set of two coupled equations: 
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These equations mimic, respectively, the conservation of 
mass and momentum for a compressible fluid of density p 
{— |Ap) and velocity V = \/(j) ; an effective stress tensor 
is given by 
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4 dxjdxh 
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where P = —^ J pf'{p)dp - is the standard hydrostatic 
pressure while the second term on the right hand side 
of Eq.(l5]) represents the quantum stress tensor. Thus, 
the Madelung fluid obeys an Euler-like equations of mo- 
tion which, due to the presence of the quantum stress 
term, is drastically different from the equations describ- 
ing the dynamics of an ideal fluid. However, in certain 
cases when, as for instance, for slowly varying perturba- 
tions, the quantum stress can be ignored and even though 
the obtained ideal fluid Euler equations are not entirely 
physically realistic, it is found to have some important 
qualitative features that are helpful in understanding the 
complex dynamics of the problems described by NSE. 

Equations ([T])-([5]) are valid for an arbitrary nonlinear- 
ity function /. For our particular case of / — p/{l+ p)'^, 
the equivalent hydrostatic pressure (P(0) = 0) may be 
written as 



P(p)=ln(l + p)^'^- 
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FIG. 1: Dynamics of initial vortices with Aoo = 2 (P > 0) 
for different t = , i = 50 and t = 300 time-moments. Left 
column corresponds to m = 1 and right column - to m = 2 
cases, respectively. Spacial patterns are presented inside out 
for clear view. 



FIG. 2: Dynamics of initial vortices with Aoo = 1.3 (P < 0) 
for different t = , t — 15 and t = 30 time-moments. Left 
column corresponds to tti = 1 and right column - to tti = 2 
cases, respectively. 



Notice that, for this pressure, the sound speed Cg = 
^/dP/dp = v^0.5p(p-l)/(H-p)3 "exists" only if 
p > 1. It should be pointed out that the above hydro- 
static pressure is positive for high densities and becomes 
negative for p < 2.16. In classical hydrodynamics [12] 
negative pressures are often associated with cavitation, 
which involves the formation of topological defects in the 
form of bubbles. For the OVS the density p = at 
r — 0, and tends to its equilibrium value Pr-j-oo — > Po > 1 
(po — ^D- Thus beyond the vortex core, the sound can 
propagate, i.e. the background is modulationally stable. 
The radial velocity of the fluid is zero while it performs 
a differential rotation with azimuthal velocity Ve = m/r. 
If 1 < Po < 2.16 , the hydrostatic pressure becomes neg- 
ative beyond the vortex core and the negative pressure 
forces may no longer be balanced by the centrifugal and 
quantum pressure forces. Thus, we expect that for the 
certain class of perturbation the cavitation could take 
place, and since p(0) remains zero due to the topological 
constrains (strictly speaking for m ^ 1) the cavitating 
bubble expanding in radial direction, might be formed. 

In order to verify the phenomenon suggested by the 



analysis of the Madelung-fluid "translation" of the NSE, 
we have carried out the stability analysis by solving 
Eq.(IT]) numerically. The dynamics of initial OVS-like 
singly- and doubly-charged field distribution was followed 
for negative P{po) < (1 < A"^ < 2.16 ) as well as pos- 
itive pressures P{po) < {A^^ > 2.16 ). For positive 
pressure (Fig.l), the initial field distributions evolve to- 
ward the stationary OVS with parameters determined by 
^oo (plotted spacial patterns are presented inside out for 
clear view). For negative pressure (Fig. 2.), on the other 
hand, no OVS are formed; instead, the vortex core ex- 
pands out constantly in time leading to the generation of 
ever expanding electromagnetic bubbles. 



We have shown that the delocalized vortex solitons as- 
sociated with an NSE, containing the particular form of 
saturating nonlincarity arising naturally in a relativistic 
pair plasma with temperature " asymmetry" , can become 
unstable for moderately intense background fields leading 
to the formation of an ever expanding core. Though the 
changes in the core plasma density and temperature are 
small, the electromagnetic fields are expelled from within 



the expanding core. The formation of such electromag- 
netic bubbles could have rather interesting implications 
for structure formation both in cosmic/astrophysical set- 
tings, and for pair plasmas that may, soon, be available 
in laboratory conditions. 
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